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Аннотация 

Введение. Освоение полярных районов Мирового океана способствовало повышению интереса к изучению 
волновых процессов в водоемах с ледяным покровом, обусловленных действием подвижной нагрузки. В 
большинстве работ отечественных и зарубежных ученых ледяной покров рассматривался как упругая или 
вязкоупругая пластина, нагруженная прямолинейно движущейся вертикальной силой. Однако при 
моделировании воздействия транспортных средств на ледяной покров представляет интерес рассмотрение 
задач, в которых сила движется по более сложной траектории. Поэтому целью данного исследования является 
разработка метода исследования поведения ледяного покрова под действием силы, движущейся по траектории 
сложной формы, подчиняясь произвольному закону движения. 

Материалы и методы. Предложен метод решения задач о действии движущейся произвольным образом силы по 
ледяному покрову водоема конечной глубины, заполненного невязкой несжимаемой жидкостью. Ледяной покров 
рассматривался как вязкоупругая пластина, лежащая на поверхности жидкости, находящейся в состоянии 
потенциального течения. На верхнюю поверхность пластины приложена сосредоточенная сила, движущаяся по 
произвольной замкнутой траектории и являющаяся периодической по времени. На нижнюю поверхность пластины 
со стороны жидкости действует гидродинамическое давление. В силу периодичности приложенной к пластине 
нагрузки при решении дифференциальных уравнений, описывающих поведение данной системы, использовали 
интегральное преобразование по времени. Далее, применяя традиционные методы, получили формулы для 
вычисления напряжений и перемещений в пластине и компонент вектора скорости частиц жидкости. Эти формулы 
представили в виде повторного интеграла. Для вычисления интегралов использовали численные методы. 
Результаты исследования. Расчеты показали, что прогиб ледяного покрова заметно увеличивается с ростом 
скорости и касательного ускорения движения нагрузки. К росту прогиба также приводят увеличение времени 
релаксации льда и уменьшение радиуса траектории движения нагрузки. Рассчитано распределение вектора 
скорости частиц жидкости по глубине водоема. 

Обсуждение и заключения. Предложенный метод показал свою эффективность при решении задач о действии 
подвижной нагрузки на ледяной покров водоема. С его помощью исследовано влияние на прогиб ледяного 
покрова кривизны траектории движения и механических свойств льда, кинематических характеристик движения 
нагрузки. Метод моделирует воздействие транспортных средств на ледяной покров водоема. Результаты его 
применения могут быть использованы при строительстве ледовых дорог или аэродромов на льду. 


Ключевые слова: бесконечный ледяной покров, движущаяся нагрузка, произвольная замкнутая траектория, 
переменная скорость. 
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Введение. В последнее время, в связи с освоением северных территорий и районов Мирового океана, огромное 
внимание уделяется проблемам изучения поведения ледяного покрова под действием внешней нагрузки. Этим 
вопросам посвящено большое количество работ отечественных и зарубежных ученых. Воздействие на 
вязкоупругую плавающую пластину импульсной движущейся нагрузки рассматривалось в [1]. Воздействие 
подвижной нагрузки на ледяной покров, примороженный к стенкам канала, представлено в [2, 3]. Колебания 
ледяного покрова, обусловленные движущейся с постоянной скоростью нагрузкой, рассматривались в [4]. 
Поведение полубесконечного ледяного покрова при действии на него равномерно движущейся нагрузки 
изучалось в [5]. Распространения волн, возбуждаемых вдоль канала с ледяным покровом, рассмотрены в 
работе [6]. В [7] при моделировании применялись нелинейные модели. В монографии [8] приведены результаты 
исследований поверхностных волн в море с плавающим битым и сплошным ледяным покровом. Разрушение 
льда под действием движущейся нагрузки рассматривалось в [9]. В приведенных выше работах в основном 
рассматривалось прямолинейное движение нагрузки. Целью данного исследования являлась разработка метода 
решения задачи о воздействии силы, движущейся по ледяному покрову произвольным образом. Достижение 
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поставленной цели позволяет более точно моделировать воздействие на ледяной покров транспортных средств, 
движение которых происходит часто по достаточно сложным траекториям и по сложному закону. 

Постановка задачи. Рассматриваются колебания бесконечного ледяного покрова, лежащего на поверхности 
водоема конечной глубины, под действием силы, движущейся произвольно по замкнутой траектории. Ледяной 
покров моделируется тонкой вязкоупругой пластиной, механические свойства которой описываются моделью 
Кельвина-Фойгта. Водоем наполнен несжимаемой жидкостью. 

Материалы и методы. Изгиб ледяного покрова описывается дифференциальным уравнением [7]: 

(1+ тои, + ЗОРИ + КИ! + Б9, Фо = 9, (1) 
где И/(х,у,ё) — прогиб ледяного покрова; д=Ей/(12(1—и?)); Е — модуль Юнга; и — коэффициент Пуассона; й — 
толщина ледяного покрова; т, — время релаксации деформаций; Д’=0,*+20,? 9. +д,*; р„ — плотность льда; р, — 
плотность воды; с? = р.й/О, К=р,э/О, = р»/О; Р(х,у,ё) — действующая на поверхности льда нагрузка; Ф(х,у, 5,1) 
— потенциал движения жидкости. Поведение жидкости подчинено уравнению: 


ДАФ=0. 
Граничные условия на границе лед-вода при <=0 и на дне водоема <=-Н (Н— глубина водоема) имели вид: 
ОИ/=0.Ф 5=0, 0.Ф &=-Н =0. 


Предполагалось, что сила Р (х, у, #) перемещается по произвольной замкнутой траектории у произвольным 
образом. При этом полагалось, что Р=Р(5(1)), где 5— дуговая координата, отсчитываемая от некоторой 
фиксированной точки кривой 7. Если функция 5(1) является периодической с периодом Т, то вполне очевидно, 
что периодической будет и величина Р(5(!)), причем также с периодом Т. Параметрическое задание траектории 
принималось в виде: 

[ = хо(®) 
у = У@)’ 
где — время. 

Рассматривая установившийся процесс и применяя интегральное преобразование по переменным х иу, а 
также конечное интегральное преобразование по { на интервале [0; 7], получим: 


(1- Фо) ро с + ИНЬ Фо. -0=Рор, 
д2Фо - р? Фо=0, 
где \» Фо — изображения неизвестных функций УМ и Ф, р’=7+ а’; ^; а — параметры интегральных 
преобразований, соответственно, по переменным х и у; ®„=2ли/Г, п=0, 1,2, 3,.... Решая второе уравнение с 


учетом граничных условий, получим: 
Фо=-ю„Тосй(р(<=Н))/р5й(РН). 

Тогда ИЗ первого уравнения следует: 

М’, =Ро/(с?+ЬсИ(рН р) в" -р*-К+Е ,т,р*). 
Движущаяся сосредоточенная сила аппроксимировалась функцией: 

Р(ху)==”ехр(-=”((х-хь (и) +(у-у,()) Ид, 
где = — параметр. 

Далее, выполняя вполне очевидные преобразования, можно получить: 


О О: п„со8 (в, (#—т)) „эт (о, (#—т)) а 
вт Г р -+К ар. я + — я 7, (рК)е р атар 


п п 


00 
"и=(с?+Ьсй (РНУ/р) в?-р*-К , г=оитор*, В? =6?+ В, 
б=х,(т)-х, В= уо(т)У. 

Зная М, можно по известным соотношениям определить формулы для вычисления компонентов вектора 
перемещения и тензора напряжений в любой точке ледяного покрова. 

При численной реализации предлагаемого метода возникает проблема, связанная с необходимостью 
вычисления интегралов, один из которых в качестве подынтегральной функции имеет сильно осциллирующую 
функцию при достаточно больших и, а второй является несобственным с бесконечным верхним пределом. 

Для вычисления интегралов от сильно осциллирующих функций применялась квадратурная формула, 
полученная на основе использования метода кубических сплайнов [11]: 


Ь М- ®х,1 ох 
Ре 1 е Н-е ' 
[его =-—, У, - (м,„-М,), 
а © 1= 7 
где А; — длины элементарных отрезков, на которые разбивается интервал [а; 5]; $(%) — аппроксимация 


Хх)кубическим сплайном, М; = $ (х). 
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со 
При вычислении несобственного интеграла использовалось приближенное соотношение т (фар = 

А о 
у ЕФ)ар при достаточно большом А. Значение А выбиралось таким, чтобы оценка допускаемой ошибки 


со < 
| ] ы Г (р)ар| была достаточно малой. 
Используя традиционные методы, можно получить необходимые оценки для вычисляемых величин. 
Например, для величины прогиба ледяного покрова данная оценка имеет вид: 


(с? +ьет( АН)! А) + А* (1+ ,л,) + 
(с? чьем (АН) А) — А*- К] +[в,л,А* ] 


А со 
\ = (фар + |, м(ф)ар. 
При вычислении суммы ряда ДЛЯ ускорения его сходимости применялся метод сигма-множителей Ланцоша. 


[ее] 


|№(2)4Ф 


А 


< 45?е^ = 


Результаты исследования. Расчеты проводились для случая действия единичной сосредоточенной силы, 
которая двигалась по замкнутой кривой, представленной на рис. 1: 


Г,:51 =. 
5(0 = ее 7» +— 
251(9 1) 2 


При этом полагалось, что толщина ледяного покрова В=0,25 м, модуль Юнга материала пластины 
Е=500 000 000 Н/м*, коэффициент Пуассона и=1/3, плотность льда р=900кг/м*, плотность воды р=1000 кг/м*, 


, “=т, 1=[0:Т]. 


глубина водоема Н=5м, &=2,5. Радиусы, определяющие форму траектории движения силы, принимались 
равными А,=15 м, А2=9 м, Ю;=3 м (рис. 1). 


У, м 


-4 -3 -2 -1 0 1 2 3 4 Х,м 
Рис. 1. Траектория движения сосредоточенной силы 


На рис. 2а приведено изменение прогиба ледяного покрова, а также максимальных значений нормальных 


напряжений (при 2=+5) 5х и 5у (рис. 6, в) и при скорости движения единичной сосредоточенной силы 
у =7,4022 м/с и касательном ускорении \=0. 


М, мх 108 2! 9х, НИМ? у, Нар 
5 | | 200 | 
200 | | | | 
0 о 2857 | 
7 -200. ге 5 
| й 
- -200 1 
ь ‚м -50 1 | 
/ Ра | У, м % 
маш и’ о 270. и. ЗВ 
й Хм. Хм -50 Х,м `-50 0 | 
а) 6) в) 


Рис. 2. Изменение прогиба И напряжений ледяного покрова: а — изменение прогиба; б — изменение напряжения 5х; 
в — изменение напряжения 5, 
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На указанных графиках красной точкой обозначено место приложения силы, а синим цветом — траектория 
движения. 

На рис. 3 представлена зависимость изменения максимального значения прогиба ледяного покрова от 
величины скорости движения силы (а) в момент {=Т/2. Положение силы на траектории в данный момент 
отмечено красной точкой (рис. 1), при этом касательное ускорение \=0. Приведено два графика, один из 
которых соответствует радиусу траектории В ‚=15 м (сплошная линия), а второй — радиусу В =5 м (пунктирная 
линия). На графике б на рис. 3 изображена зависимость максимального значения прогиба ледяного покрова от 
величины касательного ускорения в момент 1=Т в точке траектории (2;0) на рис.1 при скорости у=0. На 
графике на рис. 3 изображена зависимость максимального значения прогиба от величины времени релаксации 


деформации. 
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Рис. 3. Зависимости прогиба ледяного покрова от: а — от скорости; б — касательного ускорения; в — времени релаксации 


Поведение жидкости представлено на рис. 4, где изображено распределение вектора скоростей частиц 
жидкости по глубине водоема при {=Т/2, скорости движения силы у=3,7011 м/с и ускорении \=0. 
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Рис. 4. Распределение вектора скоростей частиц жидкости по глубине водоема 


Обсуждение и заключения. Исследовано влияние кривизны траектории движения и механических свойств 
льда на прогиб ледяного покрова скорости и ускорения движения нагрузки. Расчеты показали, что прогиб 
ледяного покрова существенно зависит от скорости движения силы и ускорения движения. 

Однако характер качественного изменения перемещений и напряжений в ледяном покрове, обусловленных 
действием движущейся силы при изменении скорости и ускорения движения, менялся слабо. 

Заметное влияние на прогиб ледяного покрова оказывают механические свойства льда, в частности, время 
релаксации. 

Полученные результаты и предложенный метод решения подобных задач могут быть использованы при 
строительстве ледовых дорог или аэродромов на льду. Кроме того, предложенный метод решения показал свою 
эффективность и может быть использован при решении других подобных задач. 


А.В. Галабурдин. Метод решения задачи о движении нагрузки по ледяному покрову водоема по сложной траектории 
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